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Abstract 

We present an analysis of non-leptonic B decays in the minimally flavour- violating 
MSSM with large tan (3. We relate the Wilson coefficients of the relevant hadronic 
scalar operators to leptonic observables, showing that the present limits on the 
B s — ► H H~ and B + — ► t + v t branching fractions exclude any visible effect in 
hadronic decays. We study the transverse helicity amplitudes of B — > VV decays, 
which exhibit an enhanced sensitivity to the scalar operators, showing that even 
though an order one modification relative to the SM is not excluded in some of 
these amplitudes, they are too small to be detected at B factories. 



* Alexander- von-Humboldt Fellow 



1 Introduction 

If new particles exist at the TeV scale, then the striking absence of evidence so far for 
their virtual effects in B or K meson mixing and decay suggests that the pattern of 
flavour- changing interactions is governed by the standard- mo del (SM) Yukawa coupling 
matrices even at the TeV scale. The minimal supersymmetric SM (MSSM) with large 
ratio tan (3 of the Higgs vacuum expectation values and no new sources of flavour violation 
in the supersymmetry-breaking Lagrangian is an example of such a minimally flavour- 
violating (MFV) theory, which nevertheless may exhibit sizeable differences from the SM 
due to Higgs exchange. The leptonic B s — > fi + n~ and B + — > t + u t decays have been 
extensively studied in this model, as well as meson mixing and B — > Dtv t [1-13]. Higgs 
exchange also generates scalar four-quark operators, which contribute to non-leptonic B 
decays. The effects of scalar operators on non-leptonic B decays have been studied in 
the MSSM (not necessarily minimally flavour-violating) and a general two Higgs doublet 
model in [14-20], mostly in connection with transverse polarization in B decays to two 
vector mesons (VV), and for specific decay modes. Some of these studies find large 
deviations from SM expectations for non-leptonic decays. 

The present work is motivated by the question whether, given the present strong 
constraints from the leptonic decays, further insight on the MFV MSSM at large tan (3 
can be derived from charmless non-leptonic B decays. To this end, extending previous 
analyses, we relate directly the Wilson coefficients of the leptonic to the relevant hadronic 
scalar operators, including charged Higgs exchange effects, and calculate the hadronic 
matrix elements in QCD factorization [21,22]. We also study observables related to the 
helicity amplitudes of B — > VV, which exhibit an enhanced sensitivity to the Higgs- 
induced scalar operators. We find that the present limit on the B s — > fi + fi~ branching 
fraction, and the observation of B + — > t + v t with a branching fraction close to the SM 
expectation, exclude any visible effects in hadronic decays, but for an academic exception: 
the positive-helicity amplitude of B — > VV modes may receive order one modifications 
relative to the SM. However, this amplitude is too small to be detected at present or 
planned B factories. 



2 Scalar four-quark operators in the MSSM with 
large tan (3 

In the SM the effective Hamiltonian for charmless B decays is 

= % E x l D) + c ^ + E W< + c *&n + + h-c, (l) 

p=u,c \ i=3 / 

where D = d or s depending on the decay mode considered, and X p D ^ = V p bV* D denotes 
a product of CKM matrix elements. The conventions for the operators Qi and the 
approximations for the short-distance coefficients Ci are given in [22] . Here we only note 
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Figure 1: Four-quark interactions mediated by neutral and charged Higgs bosons. 



that the four-quark "current-current" and "penguin" operators Q^ 2 ,Q3_io are all of the 
{V - A) X (V =F A) form. 



In the MSSM new four-quark operators are generated and the coefficients of the SM 
operators are modified. We consider the large-tan f3 scenario in a set-up, where the super- 
partner particles are somewhat heavier than the electroweak gauge bosons and the Higgs 
bosons (the "decoupling limit"), such that the leading effect is due to Higgs exchange not 
only for the neutral but also for the charged current interactions, as shown in figure 1. Of 
particular interest are the flavour- changing neutral Higgs couplings to fermions, which 
originate from a loop-induced coupling of the "wrong" Higgs field H u to the down-type 
quarks, since these couplings are enhanced by several powers of tan/3 [2,8,10]. In the 
following we use the effective couplings given in [8,10] in the decoupling limit to obtain 
the short-distance coefficients of the scalar four-quark operators from tree-level Higgs ex- 
change. The coefficients are then evolved from the electroweak scale to the bottom mass 
scale m& by the renormalization group equations. The relevant Higgs-induced terms in 
the effective Hamiltonian can be written as 



similar to (1). The "current-current" operators 

Qii = (Pibi)s+p (Djp^s-p, Q p 12 = (Pib^s+p (Djp^s-P, (p = u, c) (3) 
originate from charged Higgs exchange; the "penguin" operators 

Qiz = (Dibi) s+ p (qjq^s-p, Q\ A = (A&j)s+p (Wi)s-p, (<? = d, s, b) (4) 

from the loop-induced neutral Higgs-fermion vertices. Here i,j denote colour indices 
and (qq) s ±p = q(l ± 7s)?- The CKM factors \ { p D) in (1), (2) are now assumed to be 
composed of the effective CKM matrix elements V*f that correspond to the low-energy 
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couplings. 
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Neutral Higgs exchange: b — ► Dqq transitions 

It is straightforward to assemble the short- distance coefficients from tree-level Higgs 
exchange in terms of the effective neutral Higgs couplings given in [8,10]. Combining 
a flavour-changing and a flavour- conserving coupling, we find in the large-tan j3 limit, 
where sin/3 ~ 1, 1/ cos/3 « tan/3: 

n d 1( x_ 1 m dj rh b 6 Y y? tan 3 /3 d 

° 13 M ~ 2 (1 + e 3 tan 0) (1 + e tan /3) (1 + e, tan 0) ^ ° 14 M " °" (5) 

Here 

2 ' J ~ M 2 H0 + Ml, My W 

with c a = cos a, . . .. The e-coefficients appearing in (5) are defined in [10] and denote 
the loop-induced Higgs-fermion couplings. In the large-tan j3 MSSM products e x tan (3 
can be of order one. Just as in the b — > Dijij transitions, the coefficients of the hadronic 
Higgs penguin operators are strongly enhanced by the factor tan 3 (3. The quark masses 
fh q are the MS masses in the low-energy effective theory at the matching scale fin- 

Higgs exchange generates (Db)s-p {qq)s+p operators as well, but in this case the 
factor rhb in (5) is replaced by %, which is at most rh s . The remaining two helicity 
combinations have short-distance coefficients multiplied by a function T^j, which van- 
ishes in the large-tan/? limit. Thus, it is sufficient to consider the operators Q\ 3 14 . The 
neutral Higgs coupling to up-type quarks (second diagram in figure 1) is suppressed at 
large tan (3 relative to the down-type quarks, thus q = dj = d, s,b. In fact, the operators 
Qf 3 i4 might also be dropped due to the small down-quark mass. The operator Q\ 314: 
has the largest coefficient, but it contributes to non-leptonic decays only through loops. 
Finally, we note that the double Higgs penguin diagrams (first diagram in figure 1 with 
flavour change at both vertices) are irrelevant to non-leptonic decays due to their extra 



CKM suppression. We therefore conclude that in the MFV MSSM with large tan/3, 
only a small set of scalar penguin operators Q s { 3 14 is relevant. Of these Q S {1 is absent at 
tree-level, but it is kept for the moment, since it may be generated by renormalization 
group evolution (see discussion below). 

Charged Higgs exchange 

The operators 12 arise from the third diagram in figure 1. Once again only the 
(S + P) x (S — P) Dirac structure is dominant at large tan (3. For charmless decays we 
need only the cases Uj = uj = p = u,c, and obtain 

D m b m D tan 2 /3 D 

Mfj + (1 + e ta,np) z 

Although Ci X is enhanced only by tan 2 /3, there is no loop suppression factor ey. Due to 
the factor charged Higgs exchange is relevant in practice only for b — > s transitions. 
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Renormalization group evolution 

We first discuss the evolution of the short-distance coefficients from a typical Higgs mass 
scale, which we assume to be \in = 200 GeV, to the bottom mass scale m b = 4.2 GeV, 
when penguin diagrams are neglected. Then each pair of operators (Qn, Qu), (Q13, Qu) 
evolves independently in leading logarithmic (LL) accuracy with anomalous dimension 
matrix (in units of a s /(47r)) 

^ 2= (- 6 2 J- <8) 

With a s (m b )/a s (fiH) ~ 2.13, this results in 

C°(m b )/C°(Li H ) = CUm b )/CU» H ) « 2.20, (9) 

while C^(yu), Cf A (fi) remain zero. Using the 2-loop NDR scheme anomalous dimension 
matrix (ADM) [23] we obtain 2.35 instead and Cg{m b ) /C^ H ) = C q lA (m b ) /C q 13 (n H ) « 
0.088, but since we do not have the 1-loop correction to the initial condition of the 
scalar operators at the next-to-leading logarithmic (NLL) evolution is not fully 
consistent. In any case, we conclude that the operators Q A2 , Q\ A can be neglected to first 
approximation, since their coefficient functions are suppressed by a factor 25. 

Including penguin diagrams requires to enlarge the operator basis, since the scalar 
operators mix at the LL level into the SM penguin operators as well as their "mirror" 
copies, defined by a global exchange of left- and right chiralities of the quark fields. 
For the following discussion we neglect the electroweak penguin operators, so we deal 
with the six SM operators Qi 2 , Q3-6, their mirror copies Q/ 2 , Q'3-6, and the six scalar 
operators Q p u 12 , Q13 14. The structure of the ADM reads 

/ 76x6 06x6 06x6 \ 
06x6 76x6 06x6 
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(10) 



sc— p /sc— p sc 
\76x6 76x6 76x6/ 



where 76x6 = 76x6 * s the ADM for the SM current- current and QCD penguin operators 
(equal for the mirror operators) and 7|x6 * s a block-diagonal matrix with three identical 
2x2 blocks given by 72x2 in (8): one for Qi 112 , one for Q® 3 14 , depending on the transition, 

and one for Q\ 3 14 . The matrices 76x6° P describe the mixing of the scalar operators into 
the penguin operators. We find that Q\ x 12 and Qf 3 14 mix into the mirror penguin 
operators, while only Q\ 3 14 mixes into the SM penguins. Thus [76x~6 P ] T = (0|0|T T ) and 
[76XVJT = (r T |r r |0), where 

/n n 1 11 1 \ 

(11) 

Solving the RGE equations leaves (9) unchanged, generates the mirror QCD penguin 
operators with coefficient functions 

q D (m 6 ) « -0.71 Cf M (m b ) x [C^Qjlh) + Cf 3 (ML i = 3 . . . 6, (12) 
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and modifies the SM penguin coefficients according to C L = Cf + 5Ci with 

5Ci(m b ) « -0.71 Cf M (m b ) x C b 13 (fi H ), i = 3 . . . 6. (13) 

Since the SM penguin coefficients Cf M (mb) are small numbers, the penguin-mixing ef- 
fects are small, unless the coefficient functions of the scalar operators are of order one. 
However, due to their different chiral structure, the mirror penguin operators contribute 
differently from the standard ones to the transverse helicity amplitudes in B — ► VV 
decays as discussed below. 



3 Constraints from _B< 



and B + 



T + V T 



The natural size of the loop-induced neutral Higgs couplings e ,ey,ej is of order 0.01, 
the precise values depending on MSSM parameters. Assuming M^o = 200 GeV and 
tan /3 = 50, this allows the scalar penguin operators to have coefficients of order Cf 3 ~ 
0.01, C\ 3 ~ 0.5, which are comparable to SM penguin coefficients. However, the non- 
observation of B s — > fi + fi~ implies much stronger limits on the size of the scalar four- 
quark operator coefficient functions. 

The decay B s — > fi + fi~ proceeds via an interaction similar to the first diagram of 
figure 1 except that the lower legs are replaced by a muon pair. Since the lower vertex is 
a tree-level neutral Higgs coupling, the leptonic and hadronic decay are closely related. 
For large tan/5, a single scalar operator (Db) s p (/2/i) 5 _ p , similar in structure to Q\ 3 , 
dominates the B s — > decay amplitude. Its coefficient function is given by 



1 mbm^eYyl tan 3 (3 

2 (1 + e 3 tan 0) (1 + e tan j3) 



21 ■ 



(14) 



with 
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2,J- 



(15) 



For large tan/3, and at the level of the present experimental limit, the SM contribution 
to the decay amplitude is negligible, and the branching ratio is given by 



Br(B s 



G 2 F P B m%T Bs 
8Tr(fhb + fh s ) 2 



A 



\a 



(16) 



Comparing (5) to (14), we see that we can eliminate C MM in favour of Cf 3 in the previous 
equation and turn it into 



ejtan/3) \C*g(fj, H )\ 



2v / 27r(m fe + m s )(fi H ) m dj (n H ] 

r f 5/2 1/2 ,,(s). 



m. 



Br (5, 



1/2 



;i7) 



The present experimental upper limit on the B s — > ^i + fi branching fraction is Br(B s — > 
< 5.8 • 10~ 8 at 95% C.L. [24]. Using f Bs = 240 MeV, m s (2GeV) = 90MeV and 
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Tn,b(mb) = 4.2 GeV, and evolving both quark masses to the common scale fin = 200 GeV, 
we obtain 

(1 + e tan/3) \C^(fi H )\ < 1.4 ■ 10" 4 , (1 + e 3 tan/?) \C b 13 (fi H )\ < 7.9 ■ 10~ 3 . (18) 

When eo and/or is are negative, the coefficient functions can be larger than the values 
on the right-hand side. However, the brackets multiplying the coefficient functions enter 
the relation between the quark masses and the down-type Yukawa couplings, and hence 
(1 + e 3 tan/?) cannot become very small, if the bottom Yukawa coupling is to remain 
perturbative. We allow a factor of three enhancement of the coefficient functions to be 
conservative (that is, the brackets are required to be larger than 1/3). Including the 
factor (9) from evolution to the scale m& leads to 

\C{ 3 {m b )\ < 0.001, \C b 13 (m b ) \ < 0.05, (19) 

while Cf 3 (m b ) is a factor rhd/rh s smaller than C{ 3 {mb) and therefore negligible. Thus, the 
coefficient functions of the hadronic flavour-changing neutral Higgs penguin operators 
are constrained to be a factor of 10 smaller than the above estimates derived from 
M A o = 200 GeV and tan/5 = 50. 

The short-distance coefficient C^(/i#), arising from charged Higgs exchange, can be 
related to B + — > t + v t in a similar way. Using (7) the ratio [11,25] 

Br(ff+ -> t + u t ) MSS m _ ( m 2 B tan 2 /? \ 2 
~ Br(5+ - t+v t ) sm \ L ml + 1 + 6 tan/? J ' ^ 

is expressed in terms of C^(/j,h) as 

^=(i + g°(M m| ' 1+ ^-T^ T- < 2i > 

V m D (fx H )m b {fx H ) J 

The present average of the Babar and Belle measurements of the branching fraction is 
Bi(B + -> r+z/ r ) = (1.51 ± 0.33) • 10~ 4 [26-28]. Employing the central value \V ub \ f Bd = 
7.4 • 10~ 4 GeV and assigning a conservative 50% uncertainty to the SM prediction of the 
branching fraction, the measurement constrains R TUr to lie in the range 

0.72 < R TVt < 2.40. (22) 

Concentrating on the case D = s this implies the allowed ranges 

-0.012 < (1 + e tan/3)C 1 s 1 (^) < -0.009, 

-0.001 < (1 + e tanp)C s u (fi H ) < 0.003. (23) 

The first range corresponds to the situation, where the charged Higgs contribution is 
about twice as large as the SM one, and opposite in sign. Requiring 1 + e tan/? > 1/3 
and including the RG evolution (9) results in 

-0.08 < Cl^rrib) < -0.06, or - 0.005 < C s n {m h ) < 0.018. (24) 
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The constraint from B + — > r + u T on C{ x is not as stringent as the one from B s — > /i + yU~ on 
C( 3 , but one must remember that the charged Higgs contribution to hadronic charmless 
decays must compete with the SM tree operators rather than the penguin operators. 
In addition, since [Cf-J <C 1, the contribution (12) to the mirror penguin coefficients 
remains small. These conclusions hold a fortiori for C^, which is a factor of m^/m s 
smaller than C^. 

To conclude this section, we remark that we also performed a MSSM parameter space 
scan, calculating explicitly the loop-induced e parameters subject to the experimental 
constraints from B s — > fi + fi~, B + — > t + u t , B — > Xg'y, and AM# ds . Here we also 
included the subleading scalar operators for B s —> fi + fi~, as well as the exact expressions 
for J-~2j: ?2i an d related functions. The resulting values of the short- distance coefficients 
C^, Cf 3 are in agreement with the ranges given above. 

4 Hadronic matrix elements for B — > M1M2 

To calculate the decay amplitudes of non-leptonic, charmless B decays, we employ the 
QCD factorization (QCDF) framework [21,22]. We refer to these papers for a discussion 
of the method and to [29,30] for the definitions and notation that we adopt below. 
Let us emphasize that given the constraints on the coefficient functions, a leading-order 
treatment, where QCDF is equivalent to naive factorization [31], would suffice. However, 
it takes little additional effort to include the first-order radiative corrections. 
The matrix element of the effective Hamiltonian is written as 

(M[M' 2 \H eS \B) = Af>(MiM^ + T||B>, (25) 

p=u,c 

where 7j account for vertex, penguin and spectator-scattering terms in the QCDF for- 
mula and parameterizes the weak annihilation amplitudes. We generalize the expres- 
sion given in [29] to account for the scalar amplitudes and those from the mirror QCD 
penguin operators, such that now 

U = 5 pu [a l {M 1 M 2 ) + a^M.M,)] A{[q s u][uD}) 

+ 5 pu [a 2 (M 1 M 2 ) + af 2 (M 1 M 2 )] A{[q s D][uu\) 

+ [a p 3 (M 1 M 2 ) + a'* D (M 1 M 2 )} £ A{[q s D][qq\) 

q=u,d,s 

+ [a p A (M l M 2 ) + a' 4 pD (M 1 M 2 )} £ A{[q s q\[qD\) 

q=u,d,s 

+ c^ EW (M 1 M 2 ) ^e g A([g s D][gg]) + ^ EW (M 1 M 2 ) £ | e q A([q s q][qD}) 

q=u,d,s q=u,d,s 

+ <(MiM 2 ) A([q s D][qq}) + ^ a^M^) A([q s q][qD}) . (26) 

q=d,s q=d,s 
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The new contributions are encoded in a® ll2 (charged Higgs effects), a' p ± (mirror QCD 
penguins) and a p q4q (neutral Higgs effects), as well as modifications of the standard QCD 
penguin amplitudes a p 4 . A similar generalization applies to the annihilation amplitudes. 
Our aim is to compare the new coefficients to those present in the SM for PP, PV, VP, 
VV (P pseudoscalar, V vector meson) final states. Note that for VV, (25) and (26) 
apply separately to each of the three independent helicity amplitudes h — 0, — , +, but 
the helicity label is suppressed in our notation. 

In (26) ^4([<7mi<ZMi] [9m 2 9m 2 ]) refers to a product of decay constant, form factor and 
other factors [29,30], and the arguments indicate the flavour content of the final state 
mesons M X M 2 . Since V ± A and S ± P operators contribute differently to pseudoscalar 
and vector final states we next write^ 



a p (M x M 2 



—a p {M x M 2 ) + a' p (M x M 2 ), 


if M X M 2 


= PP, 


a' 3 p (M 1 M 2 ) + a 5 p (M 1 M 2 ), 


if M X M 2 


= PV, 


a p (M x M 2 ) - a p (M x M 2 ), 


if M X M 2 


= VP, 


-a' 3 p (M 1 M 2 )-a' 5 p (M 1 M 2 ), 


if .\ /,.!/, 


= v°v°, 


f^(a ,p (M 1 M 2 ) + a lp (M 1 M 2 )), 


if MiM 2 


= v ± v ± 



a' p (M x M 2 ) 



(M X M 2 



a' p {M x M 2 ) 



r^a' p (M x M 2 ), 

r M 2 'P 



2 , 



\MiM- 

a' p (M x M 2 )-r^a' p (M x M 2 ), 
-a' p (M x M 2 )+r^a' p (M x M 2 ), 
{ f± 1 {-a p (M x M 2 ) + rf a' p (M x M 2 )) 



if M X M 2 
if M X M 2 
if M X M 2 
if M X M 2 
if M X M 2 



PP, 
PV, 
VP, 

v°v°, 



a p q (M x M 2 ) 



Mi 



L 13q 



(M X M 2 ) 
P n q (MiM 2 ) 
X3qy M x M 2 ) 



-ff h a p {M x M, 



if M X M 2 
if M X M 2 
if M X M 2 



PP, VP, 
PV, V°V°, 



a p q (M x M 2 ) 



-a p Uq (M x M 2 ) 



a 

J± a 



iq (M x M 2 ) 
u q (M x M 2 ) 



if M X M 2 
if M X M 2 
if M X M 2 



PP, PV, 
VP, V°V°, 
V^: 



(27) 



the same relations as the last two hold between a xx and a XX £>, and between a X2 and a X2 o 



1 (0)/F±^ Ml (0) a ratio of form factors, such 



•iii 



respectively. We denote by / ± 1 

that f^ 1 ~ m B /Aq CD and f^ 1 ~ Aq CD /mB in the heavy-quark limit [30]. It follows 
that for the transverse helicity amplitudes of B — > VV decay modes the contributions 
from the new operators obey a different hierarchy in the heavy-quark limit. While in 
the SM 

A :A-:A + = 1:^:^ (28) 



m b 



T In the following we drop the superscript "_D" (referring to b — > D transitions) on the amplitude 
parameters and Wilson coefficients of the mirror penguin contributions. 
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(up to certain electromagnetic effects [32]), the Higgs contributions to the amplitude 
satisfy 



An:A-: A A 



1 : 



A 



QCD 



A 



QCD 



m 



(29) 



This effect, noted first in [16], is interesting, since it increases the sensitivity of certain 
polarization observables to the new short-distance coefficients by a factor / + 1 « 10. On 
the other hand, the absence of tensor operators implies that the formal dominance of 
the longitudinal amplitude is preserved by the Higgs contributions. 

In QCDF the a\'^ p coefficients introduced in (27) can be written at next-to-leading 
order (NLO) in the form 



a^ p (M 1 M 2 



C 



(')<f 



i±l 



Nr. 



nP(m 2 



C{±1 CpOL & 



N c 4vr 



PP P (M 2 ), (30) 



where the upper (lower) signs apply when i is odd (even). The quantities N^(M 2 ), 



V 1 { '\m 2 ), ^' ; (MiM 2 ), P l ( ' IP (M 2 ) stand, respectively, for the tree-level result ("naive 
factorization"), the 1-loop vertex correction, spectator scattering, and the penguin dia- 
grams. 

The leading-order (naive factorization) term in (30) is simply a combination of short- 
distance coefficients, except for cases where a vector meson couples to a scalar current, 
where it is zero. This is summarized by 



r(0/ 



>(')?/ 



n\'\m 2 




3,4,5,12,14s, 

6, ll,13g and M 2 = P, 

6, ll,13g and M 2 = V. 



(31) 



The NLO coefficients in (30) can mostly be expressed in terms of those already known 
from the SM operators [21,22,29,30]. For the mirror QCD penguin operators, we find 
that they are almost identical to the SM QCD penguins, that is V^(M 2 ) = Vi(M 2 ), 
H' i (MiM 2 ) = H i (M 1 M 2 ) for i = 3... 6. For the penguin contribution P^(M 2 ) one 
replaces Cj — * C[ in the SM expression and then adds the term SP^ P 6 (M 2 ) from the 
scalar operators given in (36) below. For the scalar operators, we set C® 2 and Cf 4 to 
zero (see section 2) and, using the Fierz symmetry of the NDR renormalization scheme 
for the scalar operators [23], obtain 



a xlD {M 1 M 2 ) 
a X2D {M x M 2 ) 



N r 



N r 4tt 



^ 5 (M 2 ) + ^^ 5 (M 1 M 2 ) 



(32) 
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Figure 2: Penguin contractions. Due to colour only the second diagram contributes to 
insertions of Q\ x , Q\ 3 . 



L 13q 



(M 1 M 2 ) = C q 13 N 13q , 



{M 1 M 2 



~N7 4tt 



^ 5 (M 2 ) + ^i/ 5 (M 1 M 2 ) 



(33) 



where V 5 (M 2 ), H 5 (M 1 M 2 ) can be found in [29] for PP, PV, and VP final states and in [30] 
for VV. (Although not used in the following, since we set C X2 an d Cf 4 to zero, we note that 
similarly V U (M 2 ) = V 13 (M 2 ) = V 6 {M 2 ), and H n {M 1 M 2 ) = H 13 {M 1 M 2 ) = # 6 (MiM 2 ).) 

There are no penguin contributions to (33). However, as discussed above, the in- 
sertion of scalar operators into the penguin diagrams shown in figure 2 modifies the 
evolution of the (mirror) QCD penguin operators. Accordingly, it also contributes to the 
penguin terms P%\ P (M 2 ) in a A Q (M1M2). The correction terms are proportional to the 
coefficient functions of the scalar operators and read 

,p/*, i., s G F a 



SP P (M 1 M 2 ) 

«£(MiM 2 ) 
where G* M (s) equals Gm 2 { s ) 



AttN c 
C F a s 



- C 



13 



xi 



4 m b 
- log G m \ 

3 S n Ah 



C 



13 



N 6 (M 2 )^ log ^ 
3 \i 



G 



M 2 \ 



V ± V ± with 



for M X M 2 
1 



GU S ) 



(34) 

PP,PV,VP,V Q V°, and G% h (s) [30] for 

(35) 







while G{ h (s) equals G Ma (s) for M X M 2 = PP,PV,VP, 
Similarly, for the mirror penguin coefficients 



dycj>a2(y) G(s -ie,l-y), 

and is zero for V ± V ± 



5P^(MtM 2 ) 



SP^ P (M 1 M 2 ) 



C F a s 
4ttN c 

+ C£ 

Cfc^s 
4nN c 



c 



4 w mb rJ 



(0) 



4 i mb nf ( \ 
3 lo §— - G lh( s p) 



°13 



N 6 {M 2 )l log ^-G{ h (0) 



4 
3 



N 6 (M 2 ) i log ^-G{ h (s p ) 
A 4 



(36) 



10 



where s u = 0, s c = (m c /m,b) 2 , and now G M {s) equals G^ (s) for V ± V ± and else as 

above. Note that the explicit scale dependence in 5P^\ P {M2) cancels the extra scale 
dependence of the (mirror) QCD penguin coefficients at LL accuracy. At this point we 
should mention that the constant terms in the real part of the NLO matrix elements 
should strictly speaking only be considered at the NLL order. At this order, our cal- 
culation is, however, incomplete, since we do not consider the 1-loop QCD correction 
to the initial condition of the scalar operators, and the 2-loop mixing into the penguin 
operators, as well as the small contributions from C® 2 and C\ A . Since we do not need 
precise results for the NLO terms, as will be seen below, the present approximation is 
adequate for our purpose. However, the complete NLO results for the matrix elements 
of scalar and mirror penguin operators given above might be of more general interest. 

We also calculated the weak annihilation terms originating from the scalar oper- 
ators. In some cases the annihilation amplitude can be as large as the corresponding ai 
amplitude. Since no precise estimates are needed below, we do not discuss the annihila- 
tion amplitudes further. * 

5 Non-leptonic decays 

We are now ready to discuss the question whether there are observable effects on non- 
leptonic, charmless decays due to Higgs exchange in the MSSM with large tan /3. To this 
end, we compare the new amplitudes to those present in the SM. The essential features 
can be deduced from (26). 

• Charged Higgs exchange (a^ 12 ) contributes directly to tree-dominated decays 
(such as B — > 7T7T, 7rp, pp), but must compete with the sizeable SM tree ampli- 
tudes a^2- However, since af^ 12 oc rh D , only the case of b — > suu transitions is of 
interest. But there are no tree-dominated decays of this type, since for D = s the 

fs] (s) 

tree amplitudes are doubly CKM-suppressed, A« <C Ac ■ 

• The effects from the mirror QCD penguin operators (a 3 p 4 ) must compete with 
the SM penguin amplitudes, which according to (12) requires the scalar operator 
Wilson coefficients to be of order 1 in general, and of order 0.1 in case of the 
plus-helicity amplitude in B — »■ VV . 

• The direct contribution from the FCNC Higgs couplings (a| g4? ) is an isospin- 
violating effect that must compete only with the small SM electroweak penguins, 
and is therefore most likely to lead to an observable effect. Since 0:^4 oc rh q , only 
the case q = s is of interest. For the case of D = s, the b —>■ sss transition leads 
to final states with flavour content Mi = q s s, M 2 = ss with q s the flavour of the 

*We use this occasion to point out the following corrections to [30]: The overall sign on the right- 
hand side of (A. 15) [eq.(63) in the arXiv version] must be minus. Furthermore, the expression for 
[.A3' ] in (A. 20) [eq.(68)] must contain r^ 1 — r^ 2 [r^ 1 + r^ 2 ] rather than the opposite relative sign [33]. 
(However, the unsimplificd expressions in (A. 18) [eq.(66)] arc given correctly.) 
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Table 1: Numerical results for the coefficients of some representative AS = 1 decay 
channels. The value of a\ corresponds to the SM contribution only. The final states Kr] s , 
K*fj s , Kip do not receive tree-amplitude contributions. For 0^2 and 12 , we therefore 
provide results for the final states in square brackets. 

B meson spectator antiquark. This singles out the decay modes B — > (r)('\ <p) 
and B s (i]^, <p){ri^\ <p). For the case of D = d, the potentially interesting modes 
are B — » K^K^ and B s — > K^ffi. However, in all these decays it is impossible 
to extract the EW penguin amplitude, so the new contributions must in fact be 
compared to the larger SM QCD penguins. 

We now proceed to a more detailed discussion. The numerical amplitude values given 
below depend on parameters (quark masses, form factors, etc.), for which we choose val- 
ues as given in [29,30], including some updates. Since none of our conclusions depends on 
the precise values of these parameters, we do not list them here. The Wilson coefficients 
are evaluated at the scale /i = mj = 4.2 GeV. 

B -> PP, PV 

In table 1 we show the numerical results of the a, amplitude coefficients defined in (26) 
for the decay modes B — * Krj, K*r], Kcj). (rj s in the table refers to the strange component 
of rj, see [34].) To evaluate the Higgs contributions we assume the largest values of the 
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Table 2: Numerical results for the ccj coefficients pertaining to the three helicity ampli- 
tudes of B — > VV decays. The value of a\ corresponds to the SM contribution only. To 

compare the absolute values of the helicity amplitudes the numbers for the (00, , ++) 

parameters must be multiplied by A K *^ = ^-rriBf^ (m^/l^^*, m^F^ K * , m ( j ) F^ K *). 

The estimates above use F^ Vl = 0.06 in order to compare with the maximal SM ++ 
amplitude. The final state K*<f> does not receive tree contributions. For a 12 and a\ x 12 , 
we therefore provide results for the final states in square brackets. 

coefficient functions allowed by the constraints from leptonic decays derived in section 3, 
in detail: Cf^mj,) = -0.08, Cf 3 (m b ) = 0.001, C b n {m b ) = 0.05. 

Among the Higgs penguin amplitudes a 3q is the larger of a 3qAq , since a\ q is colour- 
suppressed and is further reduced by the radiative correction given in (33). However, 
the strong constraint on C( 3 renders always negligible, in particular as it should 
be compared to the QCD penguin amplitude a\ rather than the electroweak penguin. 
This remains true for VP amplitudes despite the fact that the SM penguin amplitude is 
smaller for these final states, and for PV amplitudes, where a\ vanishes. 

B ^VV 

The effect of Higgs exchange is also negligible in case of the longitudinal amplitude in 
B — > VV decays, since it follows the same pattern as for the PV decays with M 2 = V . 
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Due to the inverted hierarchy of the transverse polarization amplitudes, see (29), the 
minus-helicity amplitude is suppressed, while the plus-helicity amplitude is enhanced 
by a factor of m^/AgcD relative to the SM. To compare the Higgs contributions to 
the plus-helicity amplitude in the SM, in table 2 we show the ctj coefficients assuming 
pB^Vi _ q Qg^ which i s the upper limit allowed in [30]. It is evident that the magnitude 
of the Higgs-induced ai coefficients is now larger for the plus amplitude than for PP, PV 
final states and the other polarization amplitudes. In fact, a\ s would now be comparable 
to the SM penguin amplitude, if it were not annihilated by the projection on the vector 
meson at tree level, see (31). Thus, among amplitudes with the same flavour topology, we 
find that only a s l2 is larger than the corresponding SM colour-suppressed tree amplitude 
«2. The mirror QCD penguin amplitude ol^ amounts to a substantial fraction of the 
standard penguin amplitude that may reach one if F+^ Vl is smaller than the assumed 
value. This would affect the azimuthal angular distribution of AS 1 = 1 decays; in practice, 
however, the effect is unobservable. Not only is the amplitude very small in absolute 
terms, but the tree amplitudes are also subleading to the penguin amplitudes in AS* = 1 
decays. 

It is straightforward to compute branching fractions, CP asymmetries and polarization 
observables including the Higgs-exchange contributions. However, since the a« param- 
eters discussed above form the basic constituents of observables, it follows that any 
modification of the SM predictions will be invisible within theoretical uncertainties. 

6 Conclusion 

Motivated by the interest in the minimally flavour- violating MSSM with large tan/? 
owing to its potentially large impact on leptonic B decays, we analyzed non-leptonic B 
decays in this model. The hadronic and leptonic flavour- changing interactions are closely 
related, which allows us to translate the present limit on the B s — > fi + fi~ branching 
fraction, and the observation of B + — > t + v t into a constraint on the Wilson coefficients 
of the relevant scalar four-quark operators. We then calculated the matrix elements 
of scalar operators and mirror QCD penguin operators at next-to-leading order in the 
framework of QCD factorization and find that the limits on leptonic B decay branching 
fractions exclude any visible effects in hadronic decays, but for an academic exception: 
the positive-helicity amplitude of B — > VV may receive order one modifications relative 
to the SM, but this amplitude is too small to be detected at present or planned B 
factories. 
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